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Abstract 

The propagation of (planar) cracks in a heterogeneous brittle material characterized 
by a random field of toughness is considered, taking into account explicitly the 
effect of the crack front roughness on the local stress intensity factor. In the so- 
called strong-pinning regime, the onset of crack propagation appears to map onto 
a second-order phase transition characterized by universal critical exponents which 
are independent of the local characteristics of the medium. Propagation over large 
distances can be described by using a simple one-dimensional description, with a 
correlation length and an effective macroscopic toughness distribution that scale 
in a non-trivial fashion with the crack front length. As an application of the above 
concepts, the arrest of indentation cracks is addressed, and the analytical expression 
for the statistical distribution of the crack radius at arrest is derived. The analysis 
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of indentation crack radii on alumina is shown to obey the predicted algebraic 
expression for the radius distribution and its dependence on the indentation load. 



Keywords: A. Crack propagation and arrest; A. Indentation and hardness; B. 
Inhomogeneous material; C. Probability and statistics; Brittle material. 



to whom correspondence should be addressed, 
e-mail: stephane . rouxSsaint-gobain . com 

2 



1 Introduction 



Brittle materials, such as ceramics and glasses, are known to be extremely sen- 
sitive to bulk and surface defects, from which cracks can be initiated eventu- 
ally leading to failure. This extreme sensitivity calls for a statistical analysis of 
crack initiation, which has been extensively developed following t he pioneering 



work of WeibuU (IWeibuU , 



19391). In this "weakest- link" approach (jFreudenthal 



19681), the analysis is focused on the initiation stage, implicitly assuming that 
the propagation stage is obtained systematically over unlimited distances, due 
to the lower loading needed for this stage as compared to the inception of 
propagation. 



However, in some cases, even though cracks are initiated, they will not propa- 
gate over large distances, and their presence may still be acceptable in service 
condition of a given structure. For example, in ceramic / metal assemblies, 
residual stresses caused by the coefficient of thermal expansion mismatch may 
prevent cracks t o traverse the bri t tle pa rt so that a weakest link hypothesis 



does not apply fjCharles and Hilc 



20021 1. In brittle-matrix composite s, crack 



arres t is also observed due to the bridging forces induced by the fibers (jEvani 



19901 ). Other examp les of such confined cracks are those produced by inden- 
tation fiLawn 



19931 ). In the latter case, the stress field has a rapid decay with 
distance from the indentation point, and hence a crack that can easily be 
nucleated may stop shortly after initiation. 



The mechanical treatment of crack arrest in a brittle material is standard: 
namely. Linear Elastic Fracture Mechanics generally holds without much re- 
striction, and hence, a crack is expected to be arrested as soon as the stress 
intensity factor K becomes smaller than the toughness K^. The purpose of 
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the present paper is to investigate the effect of a statistical distribution of 
toughness. A statistical model is considered for the steady propagation of an 
extended planar crack front through a random landscape of toughness. This 
situation is shown to give rise to a genuine phase transition at the onset 
of propagation. The recognition of this critical character opens to numer- 
ous interesting properties, not only concerning the complexity of the spatio- 
temporal development of crack propagation, but also and more importantly on 
the statistical distribution of the global stress intensity factor characterizing 
the entire front. In particular, universality instifies the independence of some 
properties on the details of the random toughness field. At a larger scale, the 
problem of crack arrest is considered in a non-uniform loading geometry. In 
this situation, some of the features revealed by the previous analysis can be 
used. As an application of the above concepts, the case of indentation cracks 
is considered since they are naturally arrested because of the decaying stress 
intensity factor. The semi-circular geometry of the crack front will allow for 
an explicit analysis of the statistical distribution of crack arrest distance. The 
results are applied on indentation of sintered alumina samples. Four different 
statistical records of indentation crack length for different loads on the same 
material are tested against the theoretical predictions. 



2 Crack pinning 



The three-dimensional problem where the crack front may be pinned at 
different locations leads to a complex problem, because elastic interac- 
tions affect the local stress intensity factor along the crack front. As a re- 
sult, the local stress intensi ty factor depend s on the ove rall front geome- 



try in a non-local manner (jBower and Ortiz , 



1991 



19931) . Curtin JCurtin, 
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1997 



19981) performed a simplified analysis of a similar problem. A mod- 



elling of planar crack propag ation in a random environ ment was also pro- 



posed by Schmittbuhl et al. (ISchmittbuhl et al. 



Fisher (Ramanathan and Fishei 



1995a|), Ramanathan and 



19981) and Bouchaud et al. (Bouchaud et al 



19931 ) that revealed a complex spatio-temporal organization of the propaga- 
tion, and scale-invariant crack front morphology. These studies show that the 
macroscopic toughness is not simply given by the average toughness and that 
tougher elements impede crack growth in a lesser extend than with a strongest 
link assumption due to the stress intensity factor enhancement caused by crack 
front curvature. One of the most striking predictions of these approaches is 
the scale invariant nature of the crack front or crack surface roughness, which 
is qualitatively (unfortunately not quantitatively) in agreement with a large 



body of experimental analysis, a review of which is presented in (jBouchauc 



19971) . 



This problem, at this level of generality, is extremely complex, and simplify- 
ing assumptions are proposed to address it. First, the crack remains planar in 
the (x, ?/)-plane where the mean crack front is aligned along the x-direction, 
whereas the y-axis is the direction of crack propagation (Fig. 1). Second, the 
distribution of local toughness is sufficiently narrow so that a first-order per- 
turbative analysis can be used to characterize the coupling of the local stress 
intensity factor K{x) to the overall front morphology h{x). Third, the elastic 
properties of the medium can be considered as homoge neous. Those assump- 



tions allow one to use a result derived by Gao and Rice (|Gao and Ricel 119891 ) 
that relates the local stress intensity factor K{x,h{x)) along the crack front 
to the crack front shape h{x), and the reference stress intensity factor, Kq, 
which would result from the same loading with a straight crack front at the 
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same mean position Hq = {h)a 



Kix, h{x)) = i^o ( 1 + - / ^^^M&^dx' ) . (1) 

71 J [X — X'Y ' 



As discussed in ( Roux et al 



20031 ). those assumptions are not yet sufficient to 
specify entirely the problem at hand. The correlation function of the toughness 
pattern comes into play to discriminate two regimes: 

• for large correlation lengths along the propagation direction, a weak-pinning 
regime is obtained, where the crack front is only mildly affected by the disor- 
der and hence the crack advance at each point along the front is continuous 
as the mean position of the crack front moves; 

• for small correlation lengths along the propagation direction, the front un- 
dergoes a series of jumps from one stable position to the next. This regime, 
referred to as "strong-pinning," is the most interesting in the sense that the 
effective toughness is expected to depart from the arithmetic average of the 
local toughness (in contrast to the weak pinning regime). 

In the following, only the latter regime is considered with a specific design of 
the local toughness distribution so that strong-pinning will be guaranteed. 



3 Propagation model 

Let us consider a planar crack front extending over the x-axis and propagating 
along the y-direction. The geometry of the front is described by h{x,t), the 
y-coordinate of the front at an abscissa x and time t. The model describes 
the propagation of the crack from one pinned configuration to the next. A 
time coordinate is introduced to label the sequence of crack front geometries. 
However, as explained below, no physical time-dependent propagation rule is 
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considered, and hence this label is used as a convenient way of characteriz- 
ing the mean crack advance. The front is discretized along the with a 
constant interval whose physical meaning can be interpreted as the correla- 
tion length of the toughness along the x-direction. The assumption that the 
crack front can be represented as a single-valued function of x is a necessary 

1). The latter impos es the 



19891). This 



assumption for the use of the perturbative result 
slope of the front to be small compared to unity ()Gao and RiceL 
in turn has consequences on the amplitude of the local toughness distribution, 
which has to be small compared to the mean value, Kc- This quantity is the 
small parameter used in all subsequent expansions. 

Under an external load that would produce a stress intensity factor -ft'o(^) for 
a straight front, Eq. (1) allows one to estimate the local stress intensity factor, 
K{x,t). The local toughness along the crack front Kc{x, h{x,t)) is assumed 
to be distributed without spatial correlation from a given distribution p{Kc). 
To make the problem non-dimensional, the scaled stress intensity factor is 
defined as A; = {K — Kc)/Kc, and the local dimensionless toughness kc = 
{Kc — Kc) I Kc- The local stress intensity factor can thus be written, as a first 
order approximation, 

k{x,t) = k^{t)^Ki{x,t), (2) 



with 



- [x — xY 



TT 



One can compute the maximum level for which the crack would not propagate, 
k^if) < kc{x, h{x,t)) — kei{x,t), for all x along the front. Thus, the onset of 
crack propagation is given by a level such that the scaled global stress intensity 
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factor assumes the value 

K{t) = min[kc{x, h{x, t)) - kei{x, t)]. (4) 

For such a level, at one point Xactit) along the front the stress intensity fac- 
tor matches exactly the local toughness. Let us now introduce the way the 
crack propagation is modeled. At the point Xact{t), the crack can overcome 
the local toughness and advance. It will directly jump to a new obstacle 
h{xacti't),t+^) = h{xact(t),t) + Sh, whcrc 5h can be either a constant length, or 
a random one. It characterizes the correlation along the propagation direction, 
and there lies the "strong-pinning" nature of the present model. At this new 
position, another local toughness is encountered, which is uncorrelated with 
the previous ones. All other sites are supposed to remain pinned during this 
elementary move of the crack front. The new front geometry and local tough- 
ness distribution is treated the same way. This means that the external loading 
is adjusted in time so that one constantly remains at the onset of propagation 
given by Eq. (4), irrespective of the variation of K{t). The physical time needed 
to achieve the elementary move of the crack front is not introduced here, and 
hence time appears here as a discrete label given to the different crack con- 
figurations. However, when inertia effects are negligible, one can reconstruct 
from the above signal K{t) what would be the evolution of the system under 
any specific history of loading in a quasi-static regime. In particular, for a 
subcritical but constant loading, /cq, starting from any configuration, one can 
compute the expected crack advance before reaching a pinned configuration 
for which nit) > ko for the first time. 

For a crack front having a finite extent (with free or periodic boundary con- 
ditions) , one can follow numerically the change of a crack front for arbitrarily 
long times, and record statistical information concerning either the crack front 
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geometry, the activity of the front with time [i.e., position of Xact{t)], or the 
loading K,{t). The latter signal is of special interest since it provides an esti- 
mate of the (reduced) critical loading n* for which an unlimited propagation 
would be obtained 

K* = max[K(t)], (5) 

and thus the latter corresponds to the macroscopic effective toughness of the 
medium, for extended propagation, while K{t) can be understood as an instan- 
taneous toughness taking into account interactions between local toughness 
values through the dependence in the crack front geometry. As shown below, 
the crack propagation can be compared to a phase transition from the pinned 
state to an unpinned {i.e., propagation) regime, where the loading plays the 
role of the control parameter. As usually observed in such second-order phase 
transition, the behavior close to k* can be characterized by a series of critical 
exponents which are universal, i.e., independent of the local disorder intro- 
duced in the rn odel, such as the distribution p{Kc) of local toughness (see 



( Sornette 



20001 ) for a presentation of the background and a review of other 
occurrences of criticality in a variety of applications including mechanical and 
geophysical examples). 



4 Statistical properties of the steady state 



After a transient regime influenced by the initial crack geometry, the model 
gives rise to a steady state during which a number of statistical observables 
can be defined and characterized. The main properties of this steady state ar e 



recalled. More details can be found in ([Tanguv et al 



1998 : 



Skoe et al. 



20021) . 
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4-1 Front roughness 



The first interesting observation concerns the crack front morphology. No spe- 
cific scale of roughness is evidenced between the discretization scale and the 
system size, but rather a scale invariant geometry appears. The crack front 
is observed to be a self-affine object such that the r.m.s. height difference 
between two points distant of Ax scales as 



{[h{x + Ax) - h{x)]YJ^ oc Ax'^, 



(6) 



with ( ^ 0.37, irres pective of the disorder. The same property can be analyzed 



with different tools dSchmittbuhl et al. 



1995b[ ). such as Fourier power spectra. 



wavelet analysis, probability density function (p.d.f.) of height differences, all 
revealing the same property. 



4-2 Activity 



The spatio-temporal distribution of activity [i.e., map of Xactif) vs. t] displays 
a highly non-trivial organization, without any characteristic scale, neither in 
space nor in time. More precisely, by studying the p.d.f. of the difference 
d = \xact{t + At) — Xactit)\ at a fixed time interval, two regimes are found: 

• at a small distance, the p.d.f. is uniform; 

• at large distances, the p.d.f. decays as a power-law d~'^ that reflects the 
kernel of the interaction term in Eq. (1). 

The cross-over distance between these two regimes, d* depends on the time 
interval At as a power-law 

d* oc A^l/^ (7) 
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where z is the so-called dynamic exponent, related to the roughness exponent 
Cby 

2; = 1 + C. (8) 
Such a behavior can be rec overed in a class of sto c hastic processes, i.e., Linear 



Fractional Stable Motion (jKrishnamurthv et al. 



20001 1 that contains both a 



wide distribution of elementary c?, and long-range temporal correlations. 



The above property reveals that the propagation consists in localized bursts of 
activity, where the clustering of activity in space and time can be decomposed 
hierarchically from the larger ones down to the smaller ones in a scale- invariant 
fashion. 



J^..3 Distribution of ejfective toughness 



The two above properties dealt with geometric organization that may be dif- 
ficult to access experimentally. However, the model also provides access to 
the loading, and hence the effective macroscopic toughness either at each in- 
stant of time K{t) or globally n*. The interesting feature is that the signal K,{t) 
contains some information about the spatio-temporal organization of activity. 

The first observation is the statistical distribution of K,{t) which shows a long 
tail of low effective toughness, but an abrupt fall-off of the distribution near 
the maximum value k*. The latter can be characterized to behave singularly 
as 

P(/€) oc (k* - nf, (9) 

irrespective of the microscopic toughness distribution p{Kc). Figure 2 shows 
the distributions -P(k) for two different distributions p{Kc), either uniform in 
[0; 1] or centered reduced Gaussian. 
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To make the connection with the activity map, it is instructive to study the 
distribution of effective toughness P{K(t),d) conditioned by a jump between 
the activity site \xact(t + 1) —Xact{t) \ = d. The clustered nature of the activity 
suggests that one needs to reach a locally strongly pinned configuration to ob- 
serve a large distance jump of the active site. Thus the distributio n P(K(t),d) 



is exp ected to concentrate progressively closer to k* as d increases fjSkoe et al. 
20021 ) ■ This is what is observed in Fig. 3. Moreover, as d becomes large com- 
pared to the discretization scale, it can be seen that the distribution P(k, d) 
converges to a unique form up to a scaling factor. Hence, one can write 

P{K,d) = d'^'^[{K* - K)d'^], (10) 

so that a characteristic correlation length can be derived for each loading 
i\\ oc (k* — k)"'^ where u = 1/a. The physical meaning of such a length scale 
is that for a subcritical loading = k, starting from any configuration in 
the steady state, the activity is localized in a region of extent The time 
needed to reach a stable configuration can be obtained from the p revious 
scaling prope r ties, At oc if< oc {k* — k)^^^ . It can be demonstrated (|Fishei . 



1985 



Kardar 



19981) that 



or this scaling property can be read off from the interaction kernel in Eq. (1). 



Knowing the statistics of jumps d distributed as a power-law of exponent —2, 
together with the scaling properties of the conditioned distribution P{K,d), 
allows one to compute the overall distribution of k, in the vicinity of k*. It 
provides the critical behavior described by Eq. (9) with 

P = u-l = j^. (12) 

Table 1 recalls the theoretical expressions of all the critical exponents in- 
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troduced so far, together with their numerical value. All exponents can be 
expressed as a function of one of them {e.g., (). However, the numerical value 
of the latter is not known theoretically. Numerical simulations are needed to 



obtain an accurate estimate. 



been used (Nattermann et al 



unctiqnal renormalization group methods have 



1993; 



Naravan and Fisher . 



1993^ to provide an 



analytic route to estimate (. However, it r emains a perturbative m ethod, which 



even when computed up to second order ()Le Doussal et al 
proximative. 



20021) , remains ap- 



5 Crack arrest 



The discussion has been focused here on the steady state propagation in a 
homogeneous environment. With those results, the problem of crack arrest 
can be addressed in a geometry where the stress intensity factor decreases 
with the propagation distance. 



Charles et al 



In fiChudnoyskv and Kunin 
I 



1987 



Jeulin ■ 



1994; 



Charles and Hild 



2002; 



2003!) ; the problem is treated by using a simple one dimensional 



description. Only the mean propagation distance {h)x is used to parameter- 
ize the crack advance. Accordingly, the stress intensity profile, K{{h)x), is 
introduced, together with some effective toughness Kc{{h)x)- At this level of 
description, the effective toughness is assumed to be given by a random func- 
tion, with a correlation length representative of the microstructure. Thus this 
random function is modeled as a piecewise constant function over a "grain 



size" 1/A, without any correlatiq r i above this scale. 



constant (jCharles and HildL 



a Poisson mosaic 



Jeulin ■ 



2002; 



Charles et al. 



'he grain size is either 



20031) or described by using 



199J). The p.d.f. for the effective toughness, /(i^c 
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is introduced. 



Let us briefly reca ll the argument developed in (jCharles and Hildl . 12002 : 



Charles et al. 



20031 ). Let F{K) = P{Kc < K) be the cumulative probability 
of effective toughness F{K) = f{Kc)dKc. The probability Q{h) that the 
crack arrest distance be larger than h is then computed by writing that the 
condition Kc{y) < K{y) hold for ally < h {i.e., a "strongest" link assumption) 



hX 



Q{h) = Y[F{K{z/X)), 



2=1 



with 



K{^/X)= mm K{y). 

j/e[(j-i)/A,«/A[ 



(13) 



(14) 



The continuous approximation of Eq. (13) 

Q{h) ^ exp (^X j\og[F{K{y))]dy^ 



(15) 



holds for long distances as compared to the "grain size" hX^ 1 and a smooth 
stress intensity factor profile XdK{y)/dy <^ K{y). This simple result was then 



exten ded to the more complex case of subcritical crack growth (ICharles et al 



20031 ). However, this situation is not considered in the present study. 



Let us note that a similar expression can be found b y assu r ning t hat the 



grain boundaries are defined by a Poisson tessellation (jjeulin 



19941 ). Under 



this assumption, 1/A is the average grain size and the propaga tion probabil 



itv with a decaying str ess intensity factor can be written as fjjeuhn 



Charles and Hild 



20021 ) 



1994 : 



Q{h) = exp (-X j\l - F{K{y))]dy\ . 



(16) 



The key question that underlies this description is the identification of the 
physical meaning of the "grain size," and the effective toughness. In this con- 
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text, the first sections that considered the roughening of the crack front are 
helpful to quantify more precisely these terms, and thus the domain of validity 
of this description. One of the major advantages of resorting to this propaga- 
tion model is that a well-defined effective toughness is now available, naturally 
parameterized in "time" t, which is proportional to the mean crack advance 
t — {h) xL / {i\ii±) where is the extension of an elementary crack advance 
during one time step {i.e., the rescaling of the length L by is used to obtain 
the number of microstructural units). 

As noted above, the crack front develops a self-affine roughness, thus the width 
of the front w = {{h^)x — (^)D^^^ dependent on the crack length L, obeying 
w oc ^_l(L/^||)^. Moreover, complex correlations in the effective toughness are 
present. However, those correlations only occur at time intervals smaller than 
the time T needed to depin the crack front over its entire length. This defines 
the larger time scale for those correlations scaling as T oc (L/^y)^. Taking into 
account the fact that L is present in the relationship between time and mean 
crack advance, the correlation length 1/A scales as 

1/A oc C±iL/C\\r-' oc C±iL/C\\)^. (17) 

Equivalently, the "grain size" 1/A is physically given by the width of the 
crack front. This signals a first non trivial feature, unexpected from the one- 
dimensional picture. The "grain size," which has to be used, depends on the 
crack front length. 

One arrives now at the most interesting output of the planar crack modeling: 
namely, the distribution of effective toughness. The power-law distribution of 
the effective toughness, which can be seen in the vicinity of k*, has already 
been discussed in Section 4.3. However, this power-law regime results from the 
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collective interaction of the crack front motion at different scales. Hence such 
a distribution holds in a domain where strong correlations are present, and 
therefore it cannot be used in the proposed approach given in Eqs. (13) and 
(15) because of the necessary assumption of statistical independence. However, 
the scale {i.e., "grain size") at which the latter condition holds is known. At 
such a scale, one has to consider the maximum effective toughness over 1/A. 
Such a quantity has already been used in the study of the effective toughness 
conditioned by the distance between consecutive active sites, d. It suffices 
to consider the largest distance d = L/2 admissible to focus on the largest 
effective toughness. Moreover, the distribution of effective toughness converges 
for large d to a universal scahng form given in Eq. (10). Thus the shape of 
the effective toughness distribution / is given by the universal function ip. 
It simply has to be translated to match k,* and scaled so that its standard 
deviation scales as {L / ^\\)~^^'^ . This is an important conclusion: When the crack 
length is large compared to the heterogeneity of the micro structure, the effective 
toughness assumes a unique universal form independent of the microstructure. 
However the standard deviation of the distribution is scale-dependent. 



Figure 4 shows the shape of the universal distribution ijj. This function has 
a power-law behavior close to the origin, with the same exponent as the one 
of the global toughness distribution, i.e., '^{x) oc x^. This implies that the 
effective toughness does not exceed the threshold k*. For large arguments 
x, 'il){x) has a rapid decay to 0. This behavior concerns the small effective 
toughness and turns out to be inessential as explicitly shown in the following 
section. 



16 



6 Application to indentation cracks 



6.1 Propagation probability 



A natural application of the above discussion concerns indentation cracks. This 
situation is a typical case where the stress intensity factor decreases strongly 
with the dist ance c from the indentation ce nter, K oc Mc~^/^, where M is the 



applied load (jPonton and Rawling; 




1989a||bO , so that cracks are generally con- 



fined to the immediate vicinity of the indent. However, radial /r nediaii crack s 



grow assuming approximately a circular or partly circular shape ()Lawn , 



19931). 



Hence the effective crack front length increases proportionally to c. This is a 
significant difference with respect to the previously treated case where L was 
assumed to remain constant. 

The change of the crack front length with c is responsible for a continuous 
change of both the length parameter 1/A and the distribution of effective 
toughness. Coming back to a one dimensional image, the problem consists 
now of a crack propagating through a series of effective grains of increasing 
size, the toughness of each grain being drawn from a distribution that is itself 
grain size dependent. The generalization of Eq. (15) leads to 



Q{c) ~ exp 







'log {FriK{r)))\{r)dr 



where the "grain size" 1/A now depends on the crack front length, 1/A = Be, 
with i? = 7r for a semi-circular crack shape, 2tt for a circular one. The B 
value is not specified in the sequel to account for any crack shape provided 
it remains self-similar at different stages of growth. Moreover, since the nu- 
merical prefactors in the scaling properties are not given, their values can be 
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disregarded at this stage 



1 . fc-^< 



and the cumulative probabihty Fr{K) of propagating through that grain de- 
rives from a distribution whose shape is universal {i.e., independent on the 
details of the toughness disorder at the microstructural scale, provided the 
latter is small compared with the crack front length) but whose width directly 
depends on the crack front length L oc c 

F,{K)^l- f{k,c)dk 

(20) 



K'-K / \ a 



dq. 



Ah 

Note that for K{c) > K*, the crack cannot be arrested so that Fc{K) — 1. This 
leads to the introduction of a characteristic crack size c* such that K[c*) = K* 
which is a lower bound on the crack size at arrest. For a large crack size as 
compared to the scale of heterogeneities, effective toughness values will always 
remain in the immediate vicinity of K* [i.e., F close to unity) so that the 
power-law behavior of -0 at the origin can be used to give a more precise 
description of \og{Fc{K)) 

log(F.(X)) oc -- [-^) . (21) 

where the scaling relations summarized in Table 1 have been used. Gathering 
these expressions, Eq. (18) can be rewritten as 

«c,.expr^(e-)-/74V"f^^^)*?^ (22, 



c- \c*J \ K* 

where A is a prefactor whose dependence with respect to the material param- 
eters ^ is 

A (X (23) 
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A change of variable x = c/c* using the scahng of the stress intensity factor 
K{c) = K*{c*/cfl'^ then leads to 

Q{xc*) ^ exp |-A(c*)2~^ u^-^ (l - _ ^24) 

The integral in Eq. (24) can be recast in terms of an incomplete beta func- 
tion. It can be noted t hat the same distribution was chosen a priori in 



( Charles and Hild 



20021 ) ■ Equation (24) therefore constitutes an a posteriori 
validation, even though the "grain size" is no longer constant in the present 
analysis. 

If the further hypothesis (m — 1) ^ 1 can be made [i.e., small standard 
deviation of the crack arrest distance as compared to c*, in practice x — 1 < 
5%), the following algebraic expression of the decay of Q{xc*) with x results 

Q{xc*) ^ exp [-A'{c*f~'^{x - 1)^^) , (25) 

where the numerical value of the above exponents is {2 — — Q ^ 2.6 

and 2 — ^ 1.63, and the dependence of A' on the characteristic pa- 
rameters and is similar to that of A [Eq. (23)]. In particular, the 
width A of the distribution of reduced crack size c/c* appears to scale as 
A oc (c*)^*-^^''-' ~ [c*)^^'^^ , becoming more and more deterministic as the 
load is increased, i.e., A oc M"^*^^"^^/^ oc M~'^''*^, with a power-law behavior 
reminiscent of a WeibuU-type result, but with the great difference that the 
exponent is independent of the concerned material! 

6.2 Analysis of indentation experiments 

The previous results are now applied to experimental data. Several indenta- 
tion tests, with different applied masses m {i.e., M oc m), have been per- 
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formed on an alumina (AI2O3) ceramic. Figure 5 shows an SEM picture of 
the material. One can see a fine-grain {i.e., lO/im) alumina polycrystal with 
an intergranular glassy phase. A microanalysis shows that the latter con- 
tains Si02, CaO and AI2O3 components. Vickers indentation with four dif- 
ferent masses have been performed on previously polished sample surfaces: 
namely, m = 0.2 kg, 0.3 kg, 0.5 kg and 1 kg. For each applied mass, mean 
indentation-generated crack lengths are measured. For the chosen loading 
range, it was checked that radial/median crack system is predominant by 
SEM observations and th at c/2a > 1, where a is the size of the mean in- 



dent half-diagonal length (jPonton and Rawlings 



1989al jbl). For each applied 



mass, measured crack lengths Cj are associated to an experimental propaga- 
tion probability Q{ci). The crack lengths are ranked in ascending order {i.e., 
ci < C2 < . . . < Cat) and the corresponding experimental probability are eval- 
uated as Q{ci) = 1 — i/{N + 1), where N is the number of measured crack 
radii for the considered applied mass. 

Equation (24) is used in the following analysis. Only two parameters have to 
be identified, namely the scale parameter A and the characteristic radius c*. 



Only one identification proced ure is presentee 



ternative methods can be used (Charles et al. 



here, although a number of al- 



20031). Following Equation (24), 



for each load series, the best c* parameter is sought by prescribing that the 
constant A be an intrinsic material parameter {i.e., unknown but identical for 
all series) . Figure 6 shows the result of the identification. A good agreement is 
obtained. Table 2 gives the value of the parameters c* for the different applied 
masses. The value of A^/*^^"^) is found to be equal to 0.15 fim~^. 

Furthermore, a rescaling procedure can be followed to compare all the exper- 
imental data on a single curve. Because of the algebraic expression of Q, one 
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has to refer to a common conventional c*, denoted as cq, to be able to compare 
all data on a unique scale. Thus a rescaled propagation probability, Qrescaied, 
is introduced such that 

Qrescaied — Q"^ (^^) 



where 

a=( , 

A check of the proposed scaling result is that Qrescaied plotted as a function of 
c/c* should fall onto a unique master curve. This is expected to be true for any 
Co choice. To treat all experimental points on an equal footing, the selection 
of the geometrical average of c* over all series was made, cq = exp((log(c*))) 
where (...) denotes the arithmetic average. Figure 7 shows the prediction for 
the indentation experiments on alumina, and a good data collapse is observed. 

Finally, Equation (24) also predicts that c* follows a power-law of the apphed 
load m. Figure 8 shows this comparison in a log-log scale together with a 
line showing the expected power- law dependence {i.e., c* oc m^/^ through 
the definition of a stress intensity factor K* oc m/(r*)^/^). Again a good 
agreement is obtained, although no more fitting parameter was used in this 
last confrontation. 



7 Conclusions 

A model for the propagation of a planar crack in a two-dimensional random 
landscape of toughness has been introduced. The complex interplay between 
the disorder and the elastic crack does not allow one to consider the successive 
failures of individual grains as independent events. The onset of propagation 
can be described as a genuine second order phase transition characterized by 
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universal critical exponents. 

When the toughness fluctuations are large enough or when the associated 
correlation length in the direction of propagation is short enough, the advance 
of the front stops being steady. The motion becomes jerky, the front jumps 
from one stable position to another one. In these strong pinning conditions the 
effective toughness of the material cannot be estimated by a simple arithmetic 
average of the toughness of the individual grains. These successive bursts 
take place over a continuous range of scales without any characteristic length 
besides the grain size and the lateral extension L of the front. Over these 
scales the front develops universal features that do not depend on the details 
of the toughness disorder. In particular, the front roughness appears to be 
self-affine: its width scales as w{L) oc L'' with C, ~ 0.37. The motion of the 
front is thus correlated up to the system size L in the lateral direction and 
in the direction of propagation. When a front of finite extent L is considered, 
a coarse-grained approach at the scale of effective grains of size L x can 
be performed. A key result is that up to scaling factors only dependent on 
the system size, the distribution of effective toughness of these effective grains 
adopts a universal form. The latter is completely independent of the details of 
the microscopic toughness. If the lateral extent of the front remains constant 
during the propagation, the latter can thus be treated as the advance of a one- 
dimensional crack through a series of effective grains characterized by the same 
well-defined effective toughness p.d.f. This result is extended to the case of self- 
similar cracks when the lateral front size is itself a growing quantity in the 
course of propagation. The coarse-grained approach remains valid but instead 
of advancing through effective grains of constant size, the one-dimensional 
crack propagates through effective grains of growing size together with a non- 
steady toughness distribution. 
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As an application to explore the predictive power of this scaling approach, 
the statistical distribution of crack arrest distances for indentation has been 
examined. The proposed scaling approach provides an analytical expression 
for the crack arrest length distribution whose dependence with respect to the 
material characteristics is limited to only two scalar parameters (a length scale, 
and a reference toughness) . The rest of the analytical expression is universal, 
(independent of the considered material). These results have been validated 
against experimental data obtained from a poly crystalline sintered alumina. 

References 

J. -P. Bouchaud, E. Bouchaud, G. Lapasset and J. Planes, Models of fractal 

cracks, Phys. Rev. Lett. 71 (1993) 2240-2243. 
E. Bouchaud, Scaling properties of cracks J. Phys. Cond. Mat. 9 [21] (1997) 

4319-4344 

A. F. Bower and M. Ortiz, A Three-Dimensional Analysis of Crack Trapping 
and Bridging by Tough Particles, J. Mech. Phys. Solids 28 [5] (1991) 815- 
858. 

A. F. Bower and M. Ortiz, The Influence of Grain Size on the Toughness of 
Monolithic Ceramics, ASME J. Eng. Mat. Technol. 115 (1993) 228-236. 

Y. Charles and F. Hild, Crack Arrest in ceramic/steel assemblies. Int. J. Fract. 
15 [3] (2002) 251-272. 

Y. Charles, F. Hild and S. Roux, Long-Term Reliability of Ceramics: the Issue 
of Crack Arrest, ASME J. Eng. Mat. Tech. 125 (2003) 333-340. 

Y. Charles, F. Hild and S. Roux, work in progress, (2003). 

A. Chudnovsky and B. Kunin, A Probabilistic Model of Brittle Crack Forma- 
tion, J. Appl. Phys. 62 (1987) 4124-4129. 



23 



W. A. Curtin, Toughening in Disordered Brittle Materials, Phys. Rev. B 55 

(1997) 11270-11276. 

W. A. Curtin, Stochastic Damage Evolution and Failure in Fiber-reinforced 

Composites, Adv. Appl. Mech. 36 (1998) 164. 
A. G. Evans, Perspectives on the Development of High- Toughness Ceramics, 

J. Am. Ceram. Soc. 73 [2] (1990) 187-206. 
D. S. Fisher, Sliding charge-density waves as a dynamic critical phenomenon, 

Phys. Rev. B 31 (1985) 1396-1427. 

A. M. Freudenthal, Statistical Approach to Brittle Fracture, in: Fracture., H. 
Liebowitz, ed., (Academic Press, New York (USA), (1968) pp. 591-619). 

H. Gao and J. R. Rice, A First Order Perturbation Analysis of Crack Trapping 
by Arrays of Obstacles, ASME J. Appl. Mech. 56 (1989) 828-836. 

D. Jeuhn, Fracture Statistics Models and Crack Propagation in Random Me- 
dia, Appl. Mech. Rev. 47 [1] (1994) 141-150. 

M. Kardar, Nonequilibrium dynamics of interfaces and lines, Phys. Rep. 301 

(1998) 85-112. 

S. Krishnamurthy, A. Tanguy, P. Abry and S. Roux, A stochastic description 
of extremal dynamics, Europhys. Lett. 51 (2000) 1-7. 

B. R. Lawn, Fracture of Brittle Solids, (Cambridge University Press, Cam- 
bridge (UK), 1993). 

P. Le Doussal, K. J. Wiese and P. Chauve, Two-loop functional renormal- 
ization group theory of the depinning transition, Phys. Rev. B 66 (2002) 
174201-1-34. 

O. Narayan and D. S. Fisher, Threshold critical dynamics of driven interfaces 

in random media, Phys. Rev. B 48 (1993) 7030-7042. 
T. Nattermann, S. Stepanow, L.-H. Tang and H. Leschhorn, Dynamics of 

interface depinning in a disordered medium, J. Phys. II (France) 2 (1993) 



24 



1483-1488. 

C. B. Ponton and R. D. Rawlings, Vickers Indentation Fracture Toughness 
Test - Part 1 - Review of Literature and Formulation of Standardized In- 
dentation Toughness Equations, Mat. Sci. Tech. 5 (1989) 865-872. 

C. B. Ponton and R. D. Rawhngs, Vickers Indentation Fracture Toughness 
Test - Part 2 - Application and Evaluation of Standardized Indentation 
Toughness Equations, Mat. Sci. Tech. 5 (1989) 961-976. 

S. Ramanathan and D. S. Fisher, Onset of Propagation of Planar Cracks in 
Heterogeneous Media, Phys. Rev. B 58 (1998) 6026-6046. 

S. Roux, D. Vandembroucq and F. Hild, Effective toughness of heterogeneous 
brittle materials, Eur. J. Mech. A/Solids [in press] (2003). 

J. Schmittbuhl, S. Roux, J. -P. Vilotte and K. J. Mal0y, Pinning of interfacial 
crack: Effect of non-local interactions, Phys. Rev. Lett. 74 (1995) 1787-1790. 

J. Schmittbuhl, J. -P. Vilotte and S. Roux, Reliability of self-affine measure- 
ments, Phys. Rev. E 51 (1995) 131-147. 

R. Skoe, D. Vandembroucq and S. Roux, Front Propagation in Random Media: 
From Extremal to Activated Dynamics, Int. J. Modern Physics C 13 (2002) 
751-757. 

D. Sornette, Critical phenomena in Natural Sciences, (Springer Verlag, Berlin 
(Germany), 2000). 

A. Tanguy, M. Gounelle and S. Roux, From individual to collective pinning: 
Effect of long-range elastic interactions, Phys. Rev. E 58 (1998) 1577-1590. 

W. WeibuU, A Statistical Theory of the Strength of Materials, (Roy. Swed. 
Inst. Eng. Res., Report 151, 1939). 



25 



List of Tables 



1 Summary of the theoretical expressions and numerical values of the critical exponents define( 



2 Values of the parameters c* for the four different applied masses m when the parameter A is 



26 



Table 1 

Summary of the theoretical expressions and numerical values of the critical 
nents defined in the text expressed as a function of C,. 



Exponent 


Expression 


Value 


c 




0.37 


u 


1/(1-0 


1.59 


a 


(1-0 


0.63 


f3 


C/(i - 


0.59 


z 


(1 + 


1.37 



Charles et al. 
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Table 2 

Values of the parameters c* for the four different applied masses m when the pa- 
rameter A is assumed to be load-independent. 



Mass m (kg) 


c* (/um) 


0.2 


11 


0.3 


14 


0.5 


21 


1.0 


36 



Charles et al. 
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Fig. 1. Depiction of crack propagation along the y-direction through a random 
toughness field shown here as a collection of inclusions which may pin the crack 
front. 
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Fig. 2. Distribution of the effective toughness K,{t) obtained for a uniform toughness 
distribution in [0; 1] (soUd Une) and centered reduced Gaussian (dashed hue). 
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Fig. 3. Distribution of the effective toughness K{t) obtained for a uniform tough- 
ness distribution in [0; 1], and conditioned distributions for geometrically distributed 
jump distances from d = 1 to d = 64 progressively concentrating close to k*. The 
conditioned distributions are not normalized so that they represent their contribu- 
tion to the overall distribution P{k). 
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Fig. 4. Universal function il>{u) giving the distribution of the effective toughness. 
Note that the origin u = corresponds to k = k*, and increasing u correspond to 
decreasing k values. 
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Fig. 5. Alumina used in the study with an average grain size of the order of 10 /xm. 
One can note the glassy intergranular phase. 
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Fig. 6. Propagation probability Q versus crack radius c for four different applied 
masses. The symbols are experimental data and the lines are identifications when 
the parameter A is assumed to be load-independent. 
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Fig. 7. Rescaled propagation probability Qrescaied versus dimensionless crack radius 
c/c* for four different applied masses. The symbols are experimental data and the 
curve is the result of the identification with a single parameter A. Prom the present 
analysis, it is expected that all experimental points should fall on the same curve. 
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Fig. 8. Parameter c* versus applied mass m. A power law with an exponent of 2/3 
fits reasonably well the experiments (dashed line). 
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